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We propose a general formula for the melting temperature of metals in terms of electronic 
mass, electronic number, and nearest-neighbor lattice distance. We derive it from the instability 
of the transverse phonon in the solid phase, using the nearly free electron model. Including 
higher order terms of vibrations enhanced near the melting temperature, the electronic restoring 
force is reduced and the ionic one is negligible. This fact greatly brings down the melting 
temperature, bringing it close to experimental data in the range of 10 % for Cs, Cu, Au, and 
Ba. Also, this theory confirms the Lindemann criterion. 
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The continuing demand for composite metals to meet 
wider range characteristics requires a fundamental and 
more practical method. Among various physical quanti- 
ties, melting temperature is one of the most fundamental 
ones. The principal requirement of a theory of melting is 
to calculate the Gibbs free energies G s and G\ of the solid 
and liquid phases as functions of pressure P and tempera- 
ture T. The melting curve in the P-T plane is then deter- 
mined by the condition G S (P,T) = Gi(P,T), where the 
calculations of G S (P,T) and G\(P,T) are taken as two 
separate problems. Statistical aspects of melting temper- 
ature have received the greatest share of attention. Mod- 
ern computing techniques have made it possible to com- 
pare various approximate schemes with one another and 
with actual and computer experiments. 1,2 ) Most studies 
have been based on inverse power law, 3-5 ) hard-core, 6 ~ 8 ) 
or other relatively idealized and short-range forces, 9, 10 ) 
such as the Lennard- Jones force field. 11 ) We cite some of 
various attempts. 12 ~ 15 ) 

Stroud and Ashcroft studied melting phenomena in 
Na based on the electron gas model with the electron- 
ion and electron-electron interactions. 16 ) They calculated 
the free energies G s and G\ of the solid and liquid phases 
separately and determined the melting temperature by 
the condition of G s — G\. After elaborate calculations, 
they obtained the melting curve which was claimed to 
be in good agreement with the results of experiments, 
up to at least 40 kbar in Na. The calculation however is 
not suitable for other metals except alkali metals, since 
in the calculation the properties of the free electron were 
fully taken into account. 

Another approach is to determine the melting temper- 
ature by the instability of the solid phase. Lindemann in- 
vestigated the instability condition of the solid phase and 
proposed a criterion for the melting temperature that the 
Lindemann ratio 5 ~ 0.1 where S is the ratio of the mean 
square amplitude of vibration of each atom about its lat- 
tice site to the nearest-neighbor distance of the lattice 
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sites. 17 ) Using the Lindemann criterion, we can obtain 
the melting temperature of any crystal. Born calculated 
the melting temperature from the vanishing point of an 
elastic stiffness constant C44, which means that the in- 
stability of the shear vibration in the solid phase occurs 
at the melting temperature. 18 ) He confirmed the Linde- 
mann criterion in the solid phase with the Lennard- Jones 
potential. Fukuyama and Platzman calculated the trans- 
verse mode instability point using the self-consistent har- 
monic approximation (SCHA) as the onset of the super- 
heating transition. They also applied the SCHA to derive 
Lindemann criterion for the alkali metals and to esti- 
mate the melting density and temperature of a Coulomb 
solid. 19,20 ) These theories, however, may not be suitable 
for metals. Metals are the most plastic solid, and the co- 
hesive energy is mainly a function of density of packing. 
Local deviations from a strict lattice regularity are easily 
accommodated. 21 ) Actually, conduction electrons derive 
an electronic potential of long-range and oscillating char- 
acter through the adiabatic principle. 

In this paper, we construct a general formula for the 
melting temperature as the vanishing point of Ct, which 
is the velocity of the transverse phonon in metals. To 
calculate c t , we use the SCHA for lattice vibrations and 
the nearly free electron model for conduction electrons. 
Finally, the melting temperature is given as 



T m = 0.145009 x 



m*R 2 k B 



(1) 



with the Boltzmann constant fcs, the Plank constant h, 
and the nearest- neighbor lattice distance i?d- The pa- 
rameters to* and n e are the effective mass and the num- 
ber of conduction electrons per site. Except for the nu- 
meric factor 0.145009, this equation is on the order of the 
Fermi temperature h 2 /(m*R^k B )^h 2 kp/(2m*kB) with 
the Fermi momentum hkp. The numeric factor is brought 
by the procedure beyond the harmonic approximation. 
Note that eq. (1) does not include the ionic mass M. 

The melting temperatures calculated using eq. (1) 
agree well with the experimental ones for alkali and no- 
ble metals. Moreover, the melting temperatures of vari- 
ous pure and composite metals can be estimated using 
the parameters as easily accessible experimental data in 
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the literature. This theory also deduces the Lindemann 
ratios of 0.183 and 0.172 for the bec and fee lattices, 
respectively. These values are consistent with the Linde- 
mann criterion. 

We proceed to microscopically derive the formula (1) 
for the melting temperature of metals. Atoms in a metal- 
lic material vibrate about their equilibrium positions. We 
denote the equilibrium position of the atom at the site I 
as Ri and the deviation as u t . Then, the ionic potential 
Vion is written as 



(2) 



l.m 



Ri - Rm and ui m = ui - u m , where 
v(Rim + «im) is the ionic potential between the atoms 
at the sites I and m. In the adiabatic approximation, the 
one-electron Hamiltonian of the conduction electrons is 

h 2 

H cl = -—V 2 + Y^w(r-R l -u l ), (3) 

where m and w(r — Ri — u{) stand for the electron 
mass and the potential of the atom at the site I with 
the coordinate Ri + ui , respectively. The electron field 
is represented as Vv(r) = J2i a i<r^(r — Ri — ui), where 
<fi(r — Ri — ui) is the Wannier function at the site I and 
aiu is the annihilation operator of the conduction elec- 
tron at the site I with the spin a. The electronic potential 
Vei is written as 

V * = (r,J C(r)#eiVv(r)^ 

= \ t ( R lm + Ulm) ((aL'w)ci + h.c.) , (4) 

l.m, a 

where t(Ri m +u lm ) = j <f>*(r - Ri - ui)H c \4>{v- R m - 
u m )dr is the transfer integral between conduction elec- 
trons at the sites I and m. (• • • ) c i denotes the thermal 
average over electronic distributions in the ionic config- 
uration {Ri + ui}. Thus, the Hamiltonian of lattice vi- 
brations is 



dui 
~~dT 



+ V ion + Vd. 



(5) 



In what follows, and respectively denote Vi on 
and V e i with ui = for all I sites, and also their differ- 



ences are AVi on = V ion - and AV ci = V ci - V cl 



(0) 



We note that {Ri} forms the minimum configuration of 

V^+V^;d(V^+V^)/dR l =0._ 

We show how to incorporate nonlinear terms to the 
restoring force and see how the procedure actually re- 
duces the restoring force. First we calculate AVS on , which 
is expanded as 

1 \ y- 1 d"r 



l,m n ai, 



n\ dRi ai ■ ■ ■ dRic 



x u lmcti ' ' ' u im.a„: 



(6) 

where d n v/dRi ai ■ ■ ■ dRi an is a derivative with ui m = 
for all I and m sites, and the suffices a\,--- ,a n 



denote components. In the SCHA, the product 
uim ai • ■ ■ ui mOLrl is decoupled in pairs. For example, the 
terms of the same type as u 2 mai u 2 ma2 u 2 ma3 reduce 

to (t)(t) u lm ai ( u lma 2 )( u lm a3 )^ whcrc (' ") denotes the 
thermal average over ionic configurations in the equilib- 
rium and the prefactor is the number of combinations 
for decoupling. Then, we apply the Fourier transforma- 
tion v(Ri) = (l/N L )J2 q v{q)e lq - Rl with N L , which is 
the number of lattice sites. Thus, AVi on is written as 



AVi on =— |p v ^ ' Ulm Y 



l,m 



x e ii<Ri-Rm) exp 



(7) 



The nonlinear terms have been incorporated into the 
exponential factor exp[— ((q ■ ui m ) 2 }/2] that reduces 
AVion- Here, we introduce a common reduction factor 
cxp[— ((q-Md) 2 )/2], where ltd is the displacement between 
the nearest-neighbor sites, since the nearest-neighbor 
terms are dominant in the equation. Defining the effec- 
tive interaction v(q) = v(q) exp[— ((q • Md) 2 )/2], we ob- 
tain 

AMon = -^J2v( q )^(iq-u lm ) 2 e i ^- R -\ (8) 

q l.m 

We apply the Fourier transformation to the displacement 
as ui = J2 q u q e lq Rl . The coefficient u q is written as 
u q = J2\ u q\ e q\ with the polarization vector e q \ for 
the mode A. We label the longitudinal mode by A = 1 
and one of the transverse modes by A = t. Then, eq. (8) 
is written as 



AVion = \ U qX U -q\ A x(G,q) 



(9) 



q,G,\ 



with A\(G, q) = v(G+q)((G+q)-e qX ) 2 - v(G)(G-e qX ) 2 , 
where G stands for the reciprocal lattice vector. 

Next, we calculate AV c i. The change in ionic config- 
uration from {Ri} to {Ri + it;} changes the potential 
in eq. (3) by Sw = w(r — Ri ui)—w(r — Ri). Accord- 
ingly, the transfer integral changes by St(ui m ) = t(Ri m + 
uim)-t(Rim), which is on the same order as 5w. The 
average of the electron operators in eq. (4) changes as 
(<4 CT <w) e i = (al^a-ma)^ + 0((Sw I e F ) 2 ) , where (• • •)^ ) 
is the thermal average for {Ri} and ep is the Fermi 
energy. Then, we obtain AV ci = (1/2) Y^i.m.a St ( u im) 



{(a\ a a ma ) 



(0) 



h.c.) in the first order of Sw/e F as an 



approximation. We manipulate AV C \ in the same man- 
ner as that for Vi on - Using the Fourier transform, a J(T = 

(l/VNE) E fe « fcCT e lfe - Ri , t(R t ) = (i/N L )j: k t(ky k - R <, 

and t(k) = t(k) cxp[— ((fe • itd) 2 )/2], we obtain 

AU cl = -L- "qAW-qA Y -fkBxik, q) (10) 



2N L 



q,\ 



with B x (k, q) = i(k + q)((k + q) ■ e qX ) 2 - i(k)(k ■ e qX f . 

Using eqs. (9) and (10), the Hamiltonian (5) exhibits 
a bilinear form. Hence, the energy uj q \ of the normal 
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vibration exhibits the dispersion 



Mu 2 qX = Y,MG,q) + w Y,hBx{k,q). (11 



The longitudinal velocity c q \ and transverse velocity c q t 
are determined using Lu q \ — c q \\q\ (e q \ || q) and Lo qt = 
Cqt 1 9 1 ( e qt -L q) for a small momentum q, respectively. 
Then, the longitudinal velocity is given as 

2 

N~l 



Mc 2 ql = v(0) + J2 r,i(G)«(G) + — ]T f k T ql (k)i(k) 

(12) 



with the operator T q \(K) = 1 - 
(l/2)(e qV K) 2 (e qV V K ) 2 iovK 
velocity is given as 



2(e„i • K)(e ql ■ V K ) + 
G or k. The transverse 



G#0 



(13) 



with the operator r, t (-K") = (1/2) (e, t ■ *0 2 (e q i ■ V K ) 2 
for K = G or fc. In the limit of q — »■ 0, we simply denote 
the velocities as ci and ct- The velocities depend on the 
renormalized ionic potential v and the transfer integral 
i, which include c\ and c t in reduction factors. Hence, we 
obtained a set of self-consistent equations (12) and (13) 
to determine c\ and ct as functions of temperature. As 
temperature increases, v and t decrease and then c\ and 
Ct decrease. We find that c t vanishes at a special temper- 
ature, which we identify the melting temperature. On the 
other hand, the longitudinal velocity c\ remains finite be- 
cause eq. (12) includes the ionic potential v(0), which is 
independent of temperature. Actually, the longitudinal 
phonon remains in the liquid phase. 

We show that near the melting point, the ionic po- 
tential v(G) for G becomes much smaller than the 
transfer integral t(k) for sc, bcc and fee crystals. The 
Fourier transforms of the ionic potential v(G) and the 
transfer energy f(fe) are given as 



v(G) = v(G)(e 
i(k) = t(k)(e 



-iG-U d \ 



v(G)e 



-ik-Ud\ 



) =t(k)e~i k «) 



(14) 
(15) 



Here, the mean square amplitude of the vibration be- 
tween the neighboring sites is calculated as 



(16) 



where Rd is the nearest-neighbor lattice coordinate and 
the last factor is calculated as 



(u q u_ q ) 



1 



2 / C t q 
9[ W 



_L ( 31. 

c\q \k B T 



N L M [c t q 

with the function g(x) = (e x — l) _1 + l/2. At high tem- 
peratures, eq. (16) reduces to 



<«?> = 



k B TR 2 d 
3M 



4 



(17) 



on the assumption of c\q/k B T, c^q/k^T <C 1. We com- 
pare two reduction factors in eqs. (14) and (15). We take 
the magnitude G m - m of one of the smallest reciprocal lat- 
tice vectors for the lower bound of \G\ in eq. (14) and 




T/Tn 



Fig. 1. c 2 J 

c to versus T/To- T m — 0.161121 is the melting temper- 
ature scaled by Tq. 



the Fermi wave number A;f for the upper bound of |fe| 
in eq. (15). Then, the ratio r(n c ) = G^ in /fc| measures 
the relative importance of v(G) against i(k), where n c is 
the electron number per site. For a simple cubic lattice, 
the ratio is written as r(n c ) = (87r/3n e ) 2 ^ 3 , which gives 
r(l) = 4.12 and r(2) = 2.60. For a bcc lattice, the ratio is 
written as r(n c ) — 2(47r/3n e ) 2 / 3 , which gives r(l) = 5.20 
and r(2) = 3.28. For an fee lattice, the ratio is writ- 
ten as r(n e ) = 3(27r/3n ) 2/3 , which gives r(l) = 4.91 
and r(2) = 3.09. These estimations show that v(G) be- 
comes negligible in comparison with i(k). Therefore, near 
the melting point, we reasonably drop the first term in 
eq. (13) and obtain 



MN 



X w Yh {(* • e g t) 2 (e„ • V fe ) 2 t(k)} ■ (18) 



This equation implies that the melting temperature is 
determined solely by the force constant given by the at- 
tractive potential due to conduction electrons. 

Since we have used the nearly free electron model, we 
have t(k) = h 2 k 2 J (2m*) with the effective electron mass 
m* . When the sound propagates along the x-direction 
with the polarization vector parallel to the y-direction, 
eq. (18) becomes 



V 



{2w) 3 MN L 



fk k. 



d 2 



v dkl \ 2m* 



h 2 k 2 



d 3 k (19) 



with k = R 2 i kBT/(9Mc 2 ) under the assumption of c t -C 
c\. By defining the constants T = 9h 2 n e / (10m* R^ks) 

and c t o = Kk^n'J 2 / (10 Mm*) 1 1 2 and introducing the 
scaled variables 7 = c 2 /c 2 and r = T/T , eq. (19) re- 
duces to 



7 



1-1 

7 



r/t 



(20) 



after integration. Note that this relation is universal, and 
in particular, M is absorbed in the scaled velocity ^/j. 

By introducing a variable x = T/7, eq. (20) is written 
as t = x(l — x)e~ x . The maximum r, r m , is determined 
by dr/dx = 0. We denote x and 7 for r = r m and 
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m* 


Ra(A) 


Ti th) CK1 


t' cx) (K) 


rTTTt 

rp(th) 

± m 






m 


j,(cx) 


Li (bec) 


1 


2.18 


3.023 


643.7 


453.7 


1.43 


Na (bec) 


1 


1.26 


3.659 


760.1 


371.0 


2.05 


K (bec) 


1 


1.25 


4.525 


501.0 


336.3 


1.49 


Rb (bec) 


1 


1.26 


4.837 


435.0 


312.6 


1.39 


Cs (bee) 


1 


1.43 


5.235 


326.9 


301.6 


1.08 


Cu (fee) 


1 


1.38 


2.56 


1418 


1358 


1.04 


Ag (fee) 


1 


1.00 


2.89 


1535 


1235 


1.24 


Au (fee) 


1 


1.14 


2.88 


1356 


1338 


1.01 


Ca (fee) 


2 


1.90 


3.95 


865.1 


1113 


0.777 


Sr (fee) 


2 


2.00 


4.30 


693.5 


1042 


0.666 


Ba (bee) 


2 


1.40 


4.35 


968.1 


1002 


0.966 


Al (fee) 


1 


1.48 


2.86 


1058.2 


933.5 


1.13 



Table I. Comparison of the theoretical melting temperature h ' 

by eq. (1) to the experimental one T^"' for various metals. Data 

(ex) 

for n e , va* , R^, and arc referred from rcf. 22. 

7 m , respectively. These values are given as 
x m = (3- Vb)/2 = 0.381966, 
r m = x m (l - x m )e- x ™ = 0.161121, (21) 
7m = r m /x m = 0.421819. 

We show 7 as a function of t in Fig. 1. We see that, 
as T increases, Ct decreases and jumps down to zero at 
T m = r m T . Thus, we arrive at eq. (1), where T m = 
0.145009 h 2 n c /(m*R 2 d k B ). 

In Table I, we compare ^ , which is the melting tem- 
perature value calculated from cq. (1), to Tm*', which is 
that from experiments, for various metals. The metals 
have monovalent and divalent elements with conductive 
s-electrons and with cubic symmetry. For to* , n c and i?<j , 
we rely on Kittel's textbook. 22 ) For Al, we replace n c = 3 
with n c = 1 , since the Hall coefficient shows that the car- 
rier in Al is a single hole. 22 -* Note that variations in M are 
typically represented by M(Cs)/M(Li) ~ 19.2 in alkali 
metals and by M(Au)/M(Cu) ~ 3.1 in noble metals. In 

(ex) 

contrast, the corresponding variations in T m ' in Table I 
are relatively small. Hence, the approximations used to 
derive the M- independent equation (1) are considered to 
be reasonable. 

We see that 7i th) agrees with Tm*' 1 in accuracy within 
10 % (Cs, Cu, Au, and Ba) and within 50 % (all metals 
in Table I except Na). In noble metals, the agreement is 
good, but in alkali metals, is somewhat higher than 
Tm X '. The general equation (1) gives definite values for 
the melting temperatures of metals without adjustable 
parameters. Although we find a certain numerical differ- 
ence between Tm h ' and Tm*' 1 , the values in Table I are 
considered to be theoretical bases for further develop- 
ment of the study of melting temperature. 

Now, we evaluate the Lindemann ratio in our theory. 
The mean square amplitude of vibration at high temper- 
atures is 



{ u ) = -^ {UqU _ q) = a— + 



(22) 



tain the Lindemann ratio S shown as 

r2 _ ( u a) a feTm 2 a 18x„ 



M 



With aR~ 2 



R\ (k F R d ) 2 

1.05224 (fee 
1.580 (bec) and 1.493 (fee), we obtain 



(23) 



- 0.05529 (bec) and 0.05224 (fee), and with 

18/(fc F i? d ) 2 



S = 0.183 (bec) and 0.172 (fee). 



(24) 



These values are consistent with and then confirm the 
Lindemann criterion S <~ 0.1. Our results also agree with 
6 — 0.15 and 0.126 obtained by numerical calculations 
with the soft core repulsive potential r~ n (n > 3) and 
the hard core repulsive potential, respectively. 8, 23 ' 

In summary, we use the shear instability of the solid 
phase to derive the equation for the melting temperature. 
We renormalize the ionic and electronic restoring forces, 
including higher order terms of displacement enhanced 
near the melting temperature by applying SCHA. The 
ionic force due to the renormalized ionic potential v(G) 
with G / decreases significantly and becomes negli- 
gible near the melting temperature. Hence, the velocity 
c t is solely determined by the force via conduction elec- 
trons. We obtain the melting temperature from the van- 
ishing point of velocity c t without any adjustable param- 
eters. The theoretical melting temperature Tm ^ is con- 

(cx) 

sistent with the experimental one T m . We can estimate 
the melting temperatures of various pure and composite 
metals if we have the values of to*, R d and n c . The Lin- 
demann criterion is obeyed. Finally, we discuss melting 
temperatures of the transition, lanthanide, and actinidc 
metals. Using the experimental data of i?d, we obtain 



A(n e /m*): A=0.7~1.2 (3d), 0.6~0.9 (4d), 



^(th) /T (ex) 
-L m / m 

and 0.8-0.5 (5d) for the transition metals; A=0. 7-0.9 
for the lanthanide metals; and A=0.5— 0.7 (Th, Pa, and 
Am) for the actinide metals. If n c /m* is about 1 — 2, 
T m th) /rm x) - 1. Actually, n c /m* - 2 in the lanthanide 
and actinide metals, and the transition metals may share 
similar values. 
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with a = (1/3N L ) Eq( 1 /? 2 )- By neglecting 1/c 2 , we ob- 13) 
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